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ABSTRACT 

Context. The Gaia-ESO Survey (GES) is a large public spectroscopic survey at the European Southern Observatory Very Large Telescope. 

Aims. A key aim is to provide precise radial velocities (^Vs) and projected equatorial velocities (vsini) for representative samples of Galactic 
stars, that will complement information obtained by the Gaia astrometry satellite. 

Methods. We present an analysis to empirically quantify the size and distribution of uncertainties in RV and vsini using spectra from repeated 
exposures of the same stars. 

Results. We show that the uncertainties vary as simple scaling functions of signal-to-noise ratio (S /N) and v sin i, that the uncertainties become 
larger with increasing photospheric temperature, but that the dependence on stellar gravity, metallicity and age is weak. The underlying uncertainty 
distributions have extended tails that are better represented by Student’s t-distributions than by normal distributions. 

Conclusions. Parametrised results are provided, that enable estimates of the RV precision for almost all GES measurements, and estimates of the 
vsini precision for stars in young clusters, as a function of S/N, vsini and stellar temperature. The precision of individual high S/N GES RV 
measurements is 0.22-0.26 km/s, dependent on instrumental configuration. 

Key words, stars: kinematics and dynamics - stars: open clusters and associations: general 


1. Introduction 

The Gaia-ESO survey (GES) is a large public survey programme 
carried out at the ESO Very Large Telescope (UT-2 Kueyen) 
with the ELAMES multi-object instrument (Gilmore et al. 2012; 
Randich & Gilmore 2013). The survey will obtain high- and 
intermediate-resolution spectroscopy of ~ 10^ stars, the majority 
obtained at resolving powers of R ~ 17 000 with the GIRAEEE 
spectrograph (Pasquini et al. 2002). The primary objectives are 
to cover representative samples of all Galactic stellar popula¬ 
tions, including thin and thick disc, bulge, halo, and stars in 
clusters at a range of ages and Galactocentric radii. The spec¬ 
tra contain both chemical and dynamical information for stars as 
faint as y ~ 19 and, when combined with complementary infor¬ 
mation from the Gaia satellite, will provide full 3-dimensional 
velocities and chemistry for a large and representative sample of 
stars. The GES began on 31 December 2011 and will continue 
for approximately 5 years. There are periodic internal and exter¬ 
nal data releases, and at the time of writing, data from the first 
18 months of survey operations have been analysed and released 
to the survey consortium for scientific exploitation - the “second 
internal data release”, known as iDR2. Part of the same data have 


* Based on observations collected with the ELAMES spectrograph at 
VLT/UT2 telescope (Paranal Observatory, ESO, Chile), for the Gaia- 
ESO Large Public Survey (188.B-3002). 


also been released to ESO through the second Gaia-ESO phase 
3 and will soon be available to the general community. 

The GES data products include stellar radial velocities (RV) 
and projected rotation velocities (vsini). A thorough under¬ 
standing of the uncertainties in RV and v sin i is an essential com¬ 
ponent of many aspects of the GES programme. Eor instance, 
the GES data are capable of resolving the kinematics of clus¬ 
ters and star forming regions, but because the RV uncertainties 
are not negligible compared with the observed kinematic dis¬ 
persion, an accurate deconvolution to establish intrinsic cluster 
velocity profiles, mass-dependent kinematic signatures, net rota¬ 
tion etc. relies on a detailed knowledge of the RV uncertainties 
(e.g. Cottaar, Meyer & Parker 2012; Jeffries et al. 2014; Lardo 
et al. 2015; Sacco et al. 2015). Searching for binary members of 
clusters and looking for outliers in RV space also requires an un¬ 
derstanding of the uncertainty distribution in order to optimise 
search criteria and minimise false-positives. Similarly, inverting 
the projected rotation velocity distribution to a true rotation ve¬ 
locity distribution (e.g. Chandrasekhar & Miinch 1950; Dufton 
et al. 2006) or comparison of the rotation velocity distributions 
of different samples requires an understanding of how uncer¬ 
tainties in V sin i broaden the observed distribution and impose 
a lower limit to the rotation that can be resolved (Erasca et al. 
2015). 
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These examples illustrate that not only does one wish to 
know the level of uncertainty in RV and v sin i as a function of 
stellar spectral type, the spectrum signal-to-noise ratio (SIN), 
the rotation rate and possibly other variables, but it is also im¬ 
portant to understand whether the uncertainties are normally dis¬ 
tributed or perhaps have extended tails that might be better rep¬ 
resented in some other way (e.g. Cottaar et al. 2014). The proce¬ 
dures for reducing and analysing the GES spectra will be fully 
detailed in forthcoming data release papers, but ultimately the 
RVs and v sin i are estimated with a detailed chi-squared fitting 
procedure (Koposov et al. in prep and Sect. 2.3). Fitting uncer¬ 
tainties can of course be computed, but these are often minor 
contributors to the overall repeatability of the measurements and 
therefore underestimate the total uncertainty. In this paper we 
empirically determine the uncertainties and their probability dis¬ 
tribution based upon repeated measurements of the same stars in 
GES. Our analysis is limited to the > 90 per cent of spectra mea¬ 
sured with the GIRAFFE spectrograph and deals only with the 
precision of the measurements, rather than their absolute accu¬ 
racy. 

In Sect. 2 we describe the GES data and the database of re¬ 
peat measurements for RV and v sin i that is available for char¬ 
acterising their uncertainties. In Sect. 3 we show how the differ¬ 
ences in RV and v sin i measured between repeated observations 
can be used to determine the underlying distribution of measure¬ 
ment uncertainty, represented by simple scaling functions that 
depend on S /N and v sin i. In Sect. 4 we investigate how these 
scaling functions alter with stellar properties. Sect. 5 considers 
how the measurement uncertainties change for different obser¬ 
vational configurations within GES. In Sect. 6, we conclude and 
provide parametric formulae and coefficients that allow an esti¬ 
mation of the RV and v sin i precision of GES measurements. 


2. Repeat measurements of radial velocities and 
projected rotation velocities 

2.1. GES observations 

The GES employs the FLAMES fibre-fed, multi-object instru¬ 
ment, feeding both the UVES high-resolution (R ~ 45 000) and 
GIRAFFE intermediate resolution (R ~ 17 000) spectrographs. 
More than 90 per cent of the spectra are obtained with GIRAFFE 
and we deal only with these data here. The Medusa fibre system 
allows the simultaneous recording of spectra from ^100 stars in 
each pointing. The stars in a single pointing are usually related 
by scientific interest (a cluster or a bulge field etc.) and cover 
a limited range of brightness (usually less than a 4 magnitude 
spread). A further ^ 15 fibres are normally allocated to patches 
of blank sky. 

The GIRAFFE spectrograph permits the recording of a lim¬ 
ited spectral range and this is selected through the use of order 
sorting filters. Eight of these have been used in the GES (HR3, 
HR5A, HR6, HRIO, HRll, HR14A, HR15N, HR21), each of 
which records a spectrum over a fixed wavelength range, al¬ 
though just three filters (HRIO, HR15N, HR21) are used for the 
large majority of observations: 

- Most observations of targets in clusters and star forming re¬ 
gions are made using order sorting filter HR15N. The wave¬ 
length range of this filter (6444-6816A) includes both the Ha 
and lithium lines and can provide useful information on the 
effective temperature (Teff), gravity (logg), age and magnetic 
activity of the target stars (Lanzafame et al 2015). 


Table 1. Numbers of short and long term repeat GIRAFFE ob¬ 
servations of RV and v sin i used for open clusters with order 
sorting filter HR15N. 


Name 

Age 

(Myr) 

Ref 

Number of repeat observations 
RV,S/N>5 vsini, >5kms“ 
short long short long 

term term term term 

Rho Ophiuchi 

1 

1 

222 

33 

34 

2 

Chamaeleon I 

2 

2 

617 

81 

108 

22 

Gamma Velorum 

6 

3 

1719 

523 

382 

80 

IC4665 

30 

4 

448 

25 

43 

1 

NGC2264 

3 

5 

2010 

333 

717 

142 

NGC2516 

140 

6 

853 

134 

266 

36 

NGC2547 

35 

7 

1045 

515 

321 

164 

NGC6633 

600 

8 

1403 

243 

103 

14 

Field giants 

— 


112 

9 

30 

2 


Notes. References. (1) Luhman & Rieke (1999); (2) Luhman (2007); 
(3) Jeffries et al. (2009); (4) Manzi et al. (2008); (5) Naylor (2009); (6) 
Meynet, Mermilliod & Maeder (1993); (7) Jeffries & Oliveira (2005); 
(8) Strobel (1991) 


- Most targets in the halo, bulge and disc fields are observed 
using both filters HRIO and HR21. The main goals here are 
to provide accurate stellar parameters and chemical abun¬ 
dances. 

GES fields are usually observed in observation blocks (OBs) 
comprising two science exposures of equal duration. In addi¬ 
tion, for filters HRIO and HR15N a short ’’simcal” exposure is 
interleaved between the science exposures. The “simcal” obser¬ 
vation illuminates five dedicated fibres with a Thorium-Argon 
(ThAr) lamp, providing a means of monitoring the wavelength 
calibration. In the HR21 observations, this role was fulfilled by 
emission lines in the sky spectra and no “simcal” exposures were 
performed. 

2.2. Data reduction 

Full details of the GES GIRAFFE data reduction will be given 
in a forthcoming paper (Lewis et al., in prep.). In brief, the raw 
data frames are corrected for a bias level using zero exposure 
bias frames and the resulting images are divided by normalised 
daytime tungsten lamp exposures to remove pixel-to-pixel sen¬ 
sitivity variations. The multiple spectra in each CCD frame are 
traced using the tungsten lamp exposures and then extracted us¬ 
ing the optimal algorithm described by Horne (1986). Given the 
readout noise and gain of the CCD, this algorithm also yields an 
estimated S /N in the extracted spectral pixels, and it is this esti¬ 
mate that is propagated through subsequent analysis steps lead¬ 
ing to the final reported S /N of the spectra. Extracted day-time 
tungsten lamp spectra are used to correct the overall shape of 
the spectrum and calibrate the individual transmission efficien¬ 
cies of each fibre. The wavelength calibration proceeded in two 
stages. Deep exposures of a daytime ThAr lamp are used to de¬ 
fine a polynomial relationship between extracted spectral pixel 
and wavelength. Then, for observations using filters HRIO or 
HR15N the wavelength calibration is modified by an offset deter¬ 
mined from the positions of prominent arc lines in the night-time 
“simcal” exposures. For observations using filter HR21 the off¬ 
set applied to the wavelength calibration is determined from the 
position of prominent emission lines in the sky spectra. Spectra 
are rebinned into 0.05 A pixels using this wavelength solution 
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and sky is subtracted using a median of the sky spectra corrected 
for the differing responses of each hbre. 

2.3. Radial velocity and projected rotation velocity estimates 

The resulting survey spectra are processed and analysed by 
working groups organised in a workflow described by Gilmore 
et al. (2012). The RV and vsini estimates used in this report 
are determined using a pipeline developed by the Cambridge 
Astronomical Survey Unit (CASU) which follows the general 
method described by Koposov et al. (2011). Details of the 
pipeline used to analyse the GES data will be described in a 
forthcoming paper (Koposov et al. in preparation). A hrst pass 
used a standard cross-correlation method with a grid of syn¬ 
thetic template spectra at a range of temperatures, metallicities 
and gravities (Munari et al. 2005) to give an initial RV estimate. 
The second pass used a direct modelling approach that hts each 
spectrum with a low-order polynomial multiplied by a template 
spectrum, with the RV, v sin i, logg, metallicity and polyno¬ 
mial coefficients as free parameters. The best fit parameter set is 
found by chi-squared minimisation with emission lines excluded 
from the htting process. The htting process is then repeated us¬ 
ing a hner grid to determine optimum values of RV and vsini 
with the other parameters held constant at their previously deter¬ 
mined values. 

The chi-squared minimisation yields an estimate of the un¬ 
certainty in the best ht parameters. However, in the case of 
GES data, this under-estimates the measurement uncertainty, in 
part due to the analysis step where spectra are re-binned but 
chiefly due to systematic uncertainties in wavelength calibration 
(Jeffries et al. 2014). Eor this reason an empirical determination 
of the measurement precision is preferred; the measurement un¬ 
certainty is estimated by comparing repeated measurements of 
RV and v sin i for the same star. 

2.4. Selected data 

To empirically characterise the RV and vsini uncertainties and 
how they depend on stellar parameters requires a database con¬ 
taining a large number of repeat observations of the same stars 
and a broad range of stellar types and rotational broadening. Eor 
these reasons, and especially to ensure a range of v sin /, we ini¬ 
tially focused on GES data for eight open clusters that were ob¬ 
served using the HR15N Alter. These clusters have ages in the 
range 1 to 600 Myr (see Table 1), covering both pre-main se¬ 
quence and main sequence objects. Only a fraction of the targets 
in each pointing will be actual cluster members, but we expect 
that cluster members will dominate any subsample of low-mass 
stars with high v sin i, since older held stars are not expected to 
rotate quickly. To provide a sample with older ages and lower 
gravities, a held consisting mainly of red giants, observed on re¬ 
peated occasions as part of the GES-CoRoT collaboration, was 
included. 

The data were restricted to observations made with two 
equal length exposures per OB. Since this is the usual mode of 
GIRAEEE observations this leads to no significant loss of data. 
Using this standard arrangement simplifies the analysis and al¬ 
lows two distinct classes of measurement uncertainty to be iden¬ 
tified: 

- Short-term repeats are where empirical estimates of uncer¬ 
tainties are obtained by comparing RV and v sin i values for 

individual targets derived from spectra measured in each of 

the individual exposures within an OB. The targets are ob¬ 


served using the same Giraffe fibre in the same configura¬ 
tion and are calibrated using the same wavelength solution. 
In this case the uncertainty is expected to be caused primar¬ 
ily by noise in the target spectra and inherent uncertainties 
in the reduction and analysis processes. Any drift in wave¬ 
length calibration over time, perhaps due to temperature or 
pressure changes, is expected to be small since the time delay 
between exposures is always <3000 s and normally <1500 s; 
there should also be no movement of the fibres and any ef¬ 
fects due to imperfect scrambling in the fibre or changing 
hour angle (see Sect. 6) should also be small. The assump¬ 
tion is also made that any significant velocity shifts due to 
binary motion on such short timescales will be rare enough 
to be neglected. 

- Long-term repeats are where uncertainties are estimated by 
comparing the mean values of RE and vsini measured in 
one OB with those measured for the same target in a sec¬ 
ond OB, where the fibre allocation and configuration on the 
plate is changed between OBs. In this case the empirical un¬ 
certainties are due to the combined effects of noise in the 
spectra, the analysis techniques plus any external uncertain¬ 
ties in the wavelength calibration or possibly differences due 
to the particular fibre used for a target or the hour angle of 
the observation. Binary motion may also contribute to any 
observed velocity shifts. A subset of these long-term repeat 
observations were observations of the same star taken on the 
same night but in a different fibre configuration. These are 
invaluable in assessing the relative importance of binaries to 
the velocity shifts. 

The data used in comparing RV measurements were se¬ 
lected to have S IN >5 (for the combined spectra in an 
OB) and those data used to compare v sin i have S /N >5 
and vsin/ > 5kms *. Table 1 shows the number of short 
and long term comparisons of RV and vsin/ available for 
each cluster. Table 2 shows the time, date, held centre co¬ 
ordinates, exposure times and numbers of targets for each 
of the Giraffe OBs used in this paper. Values of RV, v sin i, 
S /N and stellar properties are taken from the iDR2 iteration 
of analysis of the GES data, first released by the Cambridge 
Astronomical Unit to the GES working groups in May 2014 
and subsequently placed in the GES archive at the Wide 
Eield Astronomy Unit at Edinburgh UniversitjQ. 

3. Normalised distributions of measurement 
uncertainty 

Eigures 1 and 2 show the general characteristics of the ob¬ 
served RV precision, which is defined by the distribution of 
Erv - ARV/ V2, the change in RV between short-term repeat 
pairs of observations for individual targets divided by V2. Eigure 
1 shows IErvI for ~8500 short term repeats. There is a strong de¬ 
pendence on S /N and v sin i such that the measurement precision 
cannot be represented by a distribution dependent on just one of 
these parameters. Eigure 2 compares the distributions of Erv for 
short- and long-term repeats. The peak height is reduced and the 
full width half maximum (EWHM) is increased for long-term re¬ 
peats. There is thus an apparent increase in measurement uncer¬ 
tainty for targets with high S /N when compared to the precision 
assessed using short-term repeats of the same stars. 

Our general approach is to divide Erv (and the correspond¬ 
ing Eysini) by some function of the target, signal and spectro¬ 
graph properties, in order to identify the underlying normalised 

' http//ges/roe.ac.uk/ 
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Table 2. Log of VLT/Flatnes observations used in the analysis of RV and v sin i measurement precision. The full list is available as 
Supplementary Material to the on-line version of this paper. 


Filter 

Date 

observation 

UT 

RA (J2000) 
field centre 

Dec (J2000) 
field centre 

Exposure 
time (s) 

Number 

exposures 

Number 

Targets 

Cluster 

code 

HR15N 

2012-02-15 

03:07:58.00 

08:10:59.3 

-47:37:03.5 

600 

2 

111 

gam2vel 

HR15N 

2012-02-15 

03:42:56.00 

08:09:20.0 

-47:35:46.3 

600 

2 

112 

gam2vel 

HR15N 

2012-02-15 

04:18:23.00 

08:07:20.6 

-47:41:06.0 

600 

2 

81 

gam2vel 

HR15N 

2012-03-15 

03:43:29.00 

11:21:01.7 

-76:23:40.7 

600 

2 

29 

Cha-I 

HR15N 

2012-03-16 

01:39:44.00 

11:21:01.7 

-76:23:40.8 

600 

2 

29 

Cha-I 



S/N 


S/N 


Fig. 1. The empirical uncertainty in RV precision (Erv - 
ARV/ V2) estimated from the change in RV between short-term 
repeat observations of cluster targets (see Tables 1 and 2) using 
order-sorting filter HR15N. The size of the symbol indicates the 
measured value of v sin i. 



RRVH2 (km/s) 

Fig. 2. Comparison of the probability density of E^v for short- 
and long-term repeats (see Tables 1 and 2) using order-sorting 
filter HR15N. The black line shows results for short-term re¬ 
peats (i.e. pairs of observations within the same OB). The red 
histogram shows results for long-term repeats (i.e. spectra of the 
same targets but taken from different OBs where individual tar¬ 
gets are allocated to different fibres). 


distributions of measurement precision. If the underlying distri¬ 
butions are Gaussian then these normalising functions, S rv and 
Svsiai, would Correspond to the standard deviations of Erv and 
F’l.sini as a function of S/N, v sin i and stellar properties. Srv and 
S V sin i, are used here in a more general sense in order to normalise 
the Erv and Fysim distributions to an as yet unknown underlying 
distribution which could be non-Gaussian. 

Initially, we make the simplifying assumption that the nor¬ 
malising functions depend only on the S /N and v sin i of the tar¬ 
get star and on the spectrum resolution and pixel size, which are 
set by the GIRAFFE order-sorting filter. 

3.1. Normalising functions 

RV and v sin i are estimated by matching the wavelength offset 
and line width of a rotationally broadened template spectrum to 
the measured spectrum. To assess the dependency of uncertainty 
in RV on S /N and v sin i it can be shown (see Appendix A) that 
the distribution of Erv values measured from short term repeats 
scales approximately according to l{S IN) where W is the 
FWHM of individual lines in a template spectrum, rotationally 
broadened to match the line width of the measured spectrum. In 


this case (also see Appendix A), the RV precision for short term 
repeats should scale as 


S RVfl 


„(l+([vsin/]/C)2)3/4 

^- sJn - 


( 1 ) 


where C ^ 0.895c/R^, Ra is the resolving power of the spectro¬ 
graph, c is the speed of light and B is an empirically determined 
parameter that will depend on the type of star being observed. 
This is consistent with the variation of uncertainty in RV with 
S/N predicted by Butler et al. (1996) for photon limited errors. 

In the case of long-term repeats there is an additional contri¬ 
bution to the measurement uncertainty due to variations in wave¬ 
length calibration. This is independent of S /N and v sin i and 
therefore adds a fixed component A in quadrature to the short 
term uncertainty such that the distribution of Erv for long-term 
repeats scales as; 


Srv = yjA^ + Sly,, ( 2 ) 

where A will be an empirically determined constant and B and 
C are as defined in Eq. 1. 
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The relative precision of vsint used in this paper is de¬ 
fined as £„sini = Av sin ij V2{v sin i) (i.e. a fractional precision), 
where Avsint is the change between repeat observations and 
(vsint) is their mean value. To find the normalising function 
for the distribution we make the assumption that W in¬ 

creases as a function of v sin i according to the rotational broad¬ 
ening function given by Gray (1984) and that the uncertainty in 
W varies as l(S fN). In this case the uncertainty for short¬ 
term repeats (see Appendix A) scales as; 


^ „(l+([vsin/]/C)2)5/4 

(^S/N)([V Sin i]/Cf • 


(3) 


Again, a constant term is added in quadrature to account for ad¬ 
ditional sources of uncertainty present in the case of long-term 
repeats, such that the distribution of £ysini scales as; 


S 


V sin i 



2 

T sin 1,0 ’ 


(4) 


where a and f will be empirically determined constants and C is 
the same function of spectral resolution featured in Eq. 1. 


3.2. Parameters for normalising the RV measurement 
precision 

Parameters A, B and C defining the normalising function S rv are 
fitted to match the measured distribution of Erv using a dataset 
of 8,429 repeat observations, with S/N >5, taken using filter 
HR15N. Since we expect (and it turns out) that the distributions 
of these quantities are not Gaussians and have significant non- 
Gaussian tails, we choose to use the median absolute deviation 
(MAD) to characterise the observed distribution, rather than the 
square root of the mean variance which could be heavily biased 
by outliers. An estimate for the standard deviation then follows 
by noting that the MAD of a Gaussian distribution is 0.674, such 
that MAD/0.674 gives an estimate of the standard deviation. As 
we shall see, the distributions more closely follow Student’s t- 
distributions with v degrees of freedom, for which we determine 
(by Monte Carlo simulation) the corresponding corrections of 
0.82 for V = 2, 0.77 for v = 3 and 0.72 for v = 6. Uncertainties 
in the standard deviations (68 per cent confidence intervals) as a 
function of sample size are also estimated using the same Monte 
Carlo simulations. 

Defining A, B and C is then done in three steps. 

1. B is found by finding the MAD of (Erv x (S/N)/(1 + 
([v sin /]/C)^)^'^‘*, using the theoretical value of C determined 
in Appendix A (C = 15.8 kms“^ for filter HR15N, and see 
step (2) below). Figure 3a shows values of B estimated from 
data in equal bins of S /N. For S /N < 100 the average values 
per bin are close to B = 5.0km s“^ for the full data set. There 
is more scatter for S /N > 100 but the variation is not exces¬ 
sive considering the larger uncertainties due to the smaller 
numbers of data per bin. This indicates that the functional 
form of the normalising function derived in Appendix A is 
applicable to the GIRAFFE RV data. 

2. C is then checked by comparing the curve of S rv.o x (S /N), 
calculated using “empirical” values of B and C fitted to the 
measured values of Erv x (S /N) as a function of v sin i, with 
the curve predicted using B and C based on the theoretical 
value of C determined in Appendix A. Figure 3b shows that 
these two curves are very similar for the two methods, indi¬ 
cating that the theoretical value of C can be used to predict 
the scaling of measurement uncertainty in RV with v sin i . 


In fact the uncertainty on the fitted slope is largely due to the 
relatively small proportion of fast rotating stars. For this rea¬ 
son, having confirmed that the data are consistent with the 
theory in Appendix A, we prefer to use the theoretical value 
of C rather than an uncertain empirical value. The theoret¬ 
ical value for parameter C is a minimum that assumes any 
broadening of the spectral lines beyond the spectral resolu¬ 
tion is due to rotation. This is reasonable for most types of 
star in the GES, given the modest resolution of the GIRAFFE 
spectra, but if C were underestimated then we would over¬ 
estimate the increase in measurement uncertainty with v sin i 
(see Eq. 1). 

Figure 3c shows the cumulative distribution function (CDF) 
of Erv for short-term repeats normalised with S rv.o, together 
with the CDF of a unit Gaussian distribution. The distribu¬ 
tion of measurement uncertainties follows the Gaussian dis¬ 
tribution over the central region (-1 < ErvISrva < 1), but 
larger uncertainties are more frequent than predicted by the 
Gaussian. The measured distribution of ErvISrv.q is better 
represented by a Student’s t-distribution with v = 6 degrees 
of freedom, v. This value of v represents the integer value 
that provides the best fit to the normalised uncertainty of 
short-term repeats at the 5th and 95th percentiles (see Fig. 3). 
Having determined this, steps (1) and (2) are iterated, di¬ 
viding the MADs by the appropriate factor of 0.72 (for a 
Student’s t-distribution with v = 6) to estimate a true stan¬ 
dard deviation and produce the final results. 

3. The value of A that is added in quadrature to S rv.o is set to 
A - 0.25 + 0.02 km s ^ This value is chosen so that the nor¬ 
malised CDF of observed Brv found from long-term repeats, 
Erv IS rv, matches the normalised distribution of uncertainty 
from short-term repeats (ErvIS rv,o), but only between the 
upper and lower quartiles. We choose only to match this 
range because the tails of the distribution are expected to be 
different owing to the likely presence of binaries. We show in 
Sect. 4.3 that this assumption is justified because the distri¬ 
bution of ErvIS RV for those “long-term” repeats where the 
repeat observations were taken on the same observing night 
is indistinguishable from that of ErvIS rva for short-term re¬ 
peats both in the core and the tails of the distribution. 

The value of A defines the minimum level of uncertainty that 
can be achieved for GES spectra with high S IN. 

Figure 3a shows an increase in the estimated value of B for 
S IN > 100. This does not significantly affect the estimate of 
parameters A, B and C described above but does reflect the vari¬ 
ation of B with stellar properties. Fower SIN bins contain a mix 
of stars such that variations of B with stellar properties average 
out. However the smaller samples in the high S IN bins contain 
a higher fraction of stars with larger Tefr and, as we show in 
Sect. 4.1, B increases with Teff. The blue crosses in Fig. 3a show 
Teff-corrected values B as a function of S IN , using the values 
discussed in Sect. 4.1 and reported in Table 3. These show a 
more uniform variation of B with S IN. 

3.3. Parameters for normalising the vsin; precision 

Constants a, /3 and C that define the normalising function S „ sin i 
are fitted to match the measured distribution of By sin i for a subset 
of the data comprising 2004 observations with v sin i > 5 km s“'. 
Again, parameters are evaluated in three steps with the MAD 
being used to estimate the true standard deviations and the anal¬ 
ysis being iterated once the true distribution of E^smilSvsini is 
known. 
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Fig. 3. Analysis of the empirical uncertainty for short-term repeat observations of RV using filter HR15N. The solid line in plot (a) 
shows the variation of Erv x (S IN)I{1 -h ([v sin with S/N. The horizontal line indicates the value of parameter B in Eq. 1 

fitted to the full dataset. Blue crosses show the estimated values Z? as a function of S /N corrected for the measured variation of B 
with Teff (see Sect. 4.1 and Table 3). Plot (b) shows the variation of Erv x S/N with vsini. The solid line show the relationship 
predicted using the theoretical value of C and the value of B from plot (a). The dashed line shows a curve of similar functional 
form using parameters B and C htted to the binned data. In plots (a) and (b) the y-axis shows an estimate of the standard deviation 
based on the MAD divided by 0.72 (see Sect. 3.2). Plot (c) shows the cumulative probability distribution (CDF) of the normalised 
uncertainty in RV for short-term repeats. The red solid line shows results for measured data, the dashed line shows the cumulative 
distribution of a Gaussian with unit dispersion, and the diamond symbols show the cumulative distribution function for a Student’s 
t-distribution with v=6. 


1. First p is found by determining the MAD of 
£'vsini(5/A^)([vsin/]/C)2/(l -H ([vsin/]/C)2)^^'‘, using 
the theoretical value of C determined in Appendix A. The 
variation of the uncertainty with S/N shows some scatter 
(see Fig.4a) and consequently there is an +8 per cent 
uncertainty in the estimated value of (3. 

2. C is then checked by comparing the measured values of 
Ev sin i X {S /N) as a function of v sin i with the curve predicted 
using jS and C based on the theoretical value of C determined 
in Appendix A. Figure 4b shows reasonable agreement be¬ 
tween the semi-empirical curve and the measured data in¬ 
dicating that a scaling function of the form S v sin; using the 
theoretical value of C can be used to predict the variation of 
measurement uncertainty with S /N and v sin i. 

Figure 4c shows the CDF of Evsim/^vsin/.o for short-term 
repeats. This shows a more pronounced tail than the nor¬ 
malised distribution of Erv precision (Fig. 3c) such that a 
broader Student’s t-distribution with v=2 is a better ht to the 
CDF between the 5th and 95th percentiles. 

3. Finally, the value of a that represents the effect of wave¬ 
length uncertainty for long term repeats is found by match¬ 
ing the normalised £’vsini7‘5vsini distribution from 463 long¬ 
term repeats with the equivalent distribution for the short¬ 
term repeats between the upper and lower quartiles, giving 
a = 0.047 + 0.003. This corresponds to the minimum/rac- 
tional uncertainty in v sin i that can be obtained from GES 


spectra with high S /N and large v sin i. This optimum result 
is most readily achieved in spectra with vsini = 2C (i.e. 
31kms“'). Figure 4b shows that, for a given (S/N), frac¬ 
tional uncertainties rise at both higher and lower values of 
vsin/, and rise drastically for vsin/ < lOkms^^ due to the 
limited spectral resolution. 

Figure 4a shows an increase in the estimated value of /? for 
S /N > 100 due to the higher proportion of hotter stars in this 
bin. This variation is reduced when the estimated value of /S is 
corrected for the measured variation of /? with discussed in 
Sect. 4.4 and reported in Table 3. 


4. The Effect of Stellar Properties 

In Sect. 3 the constants defining the normalising functions S rv 
and S V sin / were estimated by fitting data from an inhomogeneous 
set of stars. The values obtained represent average values. In this 
section we determine how these “constants” vary with stellar 
properties, in particular T^ff, gravity, metallicity, and age. We 
make the simplifying assumption that uncertainties in RV and 
vsin / scale with S/N and vsin / as described in the last section 
and that only the parameters B in Eq. 1 and /3 in Eq. 3 depend on 
stellar properties. This follows because parameters A and a rep¬ 
resent uncertainties due to changes in wavelength calibration 
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Fig. 4. Analysis of the empirical uncertainty for short-term repeat observations of vsini using order-sorting filter HR15N. Plot 
(a) shows the variation of Fysin/ x (5/A^)([vsin/]/C)^/(l -i- ([v sinwith S/N . The horizontal line indicates the value of 
parameter (3 in Eq. 3 fitted to the full dataset. Blue crosses show the estimated values as a function of S IN corrected for the 
measured variation of p with Teff (see Sect. 4.4 and Table 3). Plot (b) shows the variation of Fysin; xSIN with v sin i. The solid line 
show the relationship predicted using the theoretical value of C and the value of p from plot (a). In plots (a) and (b) the y-axis shows 
an estimate of the standard deviation based on the MAD divided by 0.82 (for v=2, see Sect. 3.3). Plot (c) shows the cumulative 
probability distribution (CDF) of the normalised uncertainty in vsini for short-term repeats. The red solid line shows results for 
measured data, the dashed line shows the cumulative distribution of a Gaussian with unit dispersion, and the diamond symbols show 
the cumulative distribution function for a Student’s t-distribution with v=2. 


with time and fibre configuration, and parameter C should de¬ 
pend only on the spectral resolution (Eq. A4). 

4.1. Variation ofSRv with effective temperature 

Values of Tefr determined from an analysis of the iDR2 spec¬ 
tra are available in the GES archive for 75 per cent of HR15N 
targets considered in this paper. It is labelled Teffin the archive. 
The Erv values are divided between 5 evenly spaced bins of tem¬ 
perature between 3000 K and 7000 K and analysed as described 
in section 3.2. The results in Fig. 5 show a slow increase of 
B with temperature for T^ff < 5200 K such that B is within +10 
per cent of the mean value in Fig. 3b. However, above 5200 K, 
B increases rapidly with temperature to twice its mean value at 
Teff ~ 7000 K. 

The dashed lines in Fig. 5 also show results plotted as a func¬ 
tion of the ’’template” temperature (known in the GES archive as 
logTeff). This is the logarithm of the temperature of the best-fit 
synthetic spectrum that was used to determine RV and v sin i in 
the pipeline. This is likely to be less accurate than the Tefr de¬ 
rived from a full spectral analysis, but a key advantage is that 
it is available for all iDR2 targets with a RV and v sin i. In fact, 
the B values estimated using the ’’template” temperature have a 
very similar trend with and so may be used directly to es¬ 
timate temperature-dependent values of B and S rv where these 
are required. 


4.2. Variation of SRV with gravity, metaiticity and age 

Values of log g and [Fe/H] (labelled as logg and FeH in the GES 
archive) obtained from a detailed spectral analysis by the GES 
working groups are presently available for about 75 percent of 
the targets observed with order sorting filter HR15N. Analysing 
these data in bins of logg (see Fig. 6a) shows only a ~25 per¬ 
cent change in the estimated value of parameter B over a 2 dex 
range in logg. Analysis in bins of [Fe/H] (see Fig. 6b) shows 
a similarly small change in B with metallicity over the range 
of metallicities - l<[Fe/H]<l. Below this, the estimated value of 
parameter B appears to increase sharply with decreasing [Fe/H] 
but in truth there are too few data points for filter HR15N with 
[Fe/H]<-1 to estimate parameter B with any degree of accuracy. 
We confirmed that any variation seen in Fig. 6 is not due to dif¬ 
ferences in temperature - the median values of logTeff are very 
similar in all binned subsamples. 

Although the fundamental cause of any variation of RV pre¬ 
cision with age would likely be due to the evolution of logg 
in pre-main sequence stars, it is nevertheless important to con¬ 
firm that the prescription for calculating RV precision is valid 
at all ages, since studying the dynamics of young clusters is a 
key GES objective. Figure 7 shows the variation of B with stellar 
age. The adopted ages for cluster stars are those given in Table 1. 
For this plot, we attempted to separate genuine cluster members 
from field objects by selecting according to RV. For most cluster 
datasets there was a clear RV peak corresponding to the clus- 
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Fig. 5. Variation of parameter B of the scaling function for un¬ 
certainty in RV (SRv) with effective temperature. The solid line 
shows results for filter HR15N as a function of Tefr (see section 
4.1). Dashed lines show results for filters HRIO, HR15N and 
HR21 as a function of the temperature of the template spectrum 
fitted in the CASU pipeline (see Sect. 2.3). Numbers equal the 
sample size per bin. 


ter, so cluster members were selected from a range +5kms“' 
either side of this peak with little contamination. However, no 
selection by RV was made for the COROT sample or for the 
cluster NGC6633, since neither showed a clear peak in their 
RV distributions. We assume these datasets contain mostly older 
(> 1 Gyr) field stars. Figure 7 shows in any case that there is 
a weak dependence of B on age. However, it can be seen that 
this small variation is directly linked to the decreasing median 
temperatures of the cluster samples at younger ages. 

4.3. Variation ofSuv with time between observations 

In our model of RV uncertainty we assume that A represents 
some additional uncertainty arising from random changes in 
wavelength calibration with time and the effects of changes in 
fibre allocation. We fitted A, using the interquartile range of 
the uncertainty distribution in long-term repeats, in an effort 
to avoid modelling tails that might be due to binary motion. 
This simplifying assumption can be tested by plotting values of 
A determined for samples with increasing time differences be¬ 
tween observations. Figure 8a shows that the value of A depends 
only weakly on the time between observations, increasing from 
0.23 + 0.02km s“* for measurements made in different configu¬ 
rations on the same night to 0.26 + 0.02 km s ' for intervals of 
up to 100 days between observations. This confirms that the A 
value is not unduly influenced by any binaries in the sample. 
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Fig. 6. Variation of parameter B of the scaling function for un¬ 
certainty in RV (S rv) with logg and metallicity for order sorting 
filters, HRIO, HR15N and HR21. (a) shows the value of B de¬ 
termined for data in three equal bins of log g. Numbers indicate 
the sample size, (b) shows the value of B determined for data 
in three equal bins of [Fe/H]. For filter HR15N only the values 
shown in the upper two bins are reliable due to the low number 
of targets with [Fe/H]<-1. 
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Fig. 7. Variation of parameter B of the scaling function for un¬ 
certainty inRV (S rv) with target age for observations with order 
sorting filter HR15N. Square symbols show the value of B for 
targets identified as possible cluster members from their RV ver¬ 
sus the nominal age of the cluster (see Table 1). No selection by 
RV is made for NGC6633 or COROT and we assume the stars 
have mean ages > 1 Gyr. The dotted line indicates the median 
value of Teff (right hand axis values) for members identified in 
each cluster. 


The effect of binaries is far more apparent in the tails of the 
distributions. Figure 8b compares CDFs of the normalised dis¬ 
tribution of measurement uncertainty derived from the change 
in RV between short-term repeats, normalised with S rv,o (Eq- 
1), with (i) all the long-term repeats, with uncertainties nor¬ 
malised to SRV (Eq. 2); (ii) a separate distribution of ErvISrv 
for just those long-term repeats where the repeat observations 
were on the same night (nullifying the effects of all but the 
rarest, short-period binaries). By design, the three CDFs are very 
close in the interquartile range; but while the CDF for long- 
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Fig. 8. The dependence of the scaling function for uncertainty 
in RV {SRv) as a function of time between observations. Plot 
(a) shows how scaling parameter A varies with time between re¬ 
peat observations. Numbers indicate size of the sample used to 
determine A. Plot (b) shows the CDFs of the normalised distri¬ 
bution of RV precision for short-term repeats (black line), long¬ 
term repeats (red line) and long term repeats where the obser¬ 
vations were taken on the same observing night (blue dashed 
line). The CDFs for the short-term repeats and the long-term 
repeats within a night are indistinguishable using a two-tailed 
Kolmogorov-Smirnov test. 


term repeats has a more pronounced tail, better described by a 
Student’s t-distribution with v-3, the long-term repeats within 
a night are indistinguishable (with a Kolmogorov-Smirnov test) 
from the short-term repeats, following a Student’s t-distribution 
with V = 6. This is consistent with a fraction of the sample being 
binary stars that show genuine RV changes between observations 
on timescales longer than a day. It also justifies an assumption 
that the true uncertainties in a single RV measurent are best rep¬ 
resented by the v = 6 Student’s t-distribution multiplied by Srv 
as given by Eq. 2. 


Fig. 9. Variation of the scaling function for v sin i with tempera¬ 
ture and time between observations. Plot (a) shows how the fi 
parameter in Eq. 3 varies with The solid line shows re¬ 
sults using Teff from a detailed spectral analysis; the dashed 
line shows the results using the ’’template” temperature (see 
Sect. 4.4). Labels indicate the sample size per bin. Plot (b) shows 
the CDE of the normalised distributions of v sin i uncertainty for 
short and long term repeats. Also shown (as small diamonds) is 
a Student’s t-distribution with v=2. 

4.4. Variation of S^,sini with temperature and time between 
observations 

In Fig. 9a we show how (3, the parameter in the scaling function 
governing v sin i precision (see Eq.3), depends on stellar temper¬ 
ature and the time between observations . The data were divided 
into 5 equal bins of temperature. Results are shown using the 
temperature derived from detailed spectral analysis (Teff) and the 
best-fitting ’’template” temperature (logTeff). There appears to be 
little variation with below 6000 K using either temperature 
estimate, but like the parameter B governing RV precision, there 
is a rapid growth in (3 for hotter stars - by about a factor of 2 at 
Teff 7000 K. 

Figure 9b compares the CDF of the normalised measure¬ 
ment precision in v sin i for short- and long-term repeats. There 
is much less difference between these CDFs than that found be¬ 
tween the short- and long-term repeat estimates of RV precision. 
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Table 3. Constants describing the scaling function of measure¬ 
ment precision in /?y and v sin i as a function of S /N and v sin i 
(see Eqns. 2 and 4). 


Characteristics of order sorting filter 


Filter 

HRIO 

HR15N 

HR21 

Mean A (A) 

5470 

6630 

8728 

Resolution 

19800 

17000 

16200 

Range A (A) 

270 

370 

504 

Parameters defining the scaling function Srv (Eq. 3) 

A(kms ') 

0.22±0.02 

0.25±0.02 

0.26±0.02 

C(kms“') 

13.6 

15.8 

16.6 

Average value of B for the mix 

of stars analysed in this paper 

B (kms 

2.3 

5.0 

7.1 

Variation of B with template temperature (section 4.1) 

B (3200-4000 K) 

1.8±0.2 

3.9±0.1 

7.3±0.2 

B (4000-4800 K) 

1.7±0.1 

4.4±0.2 

6.8±0.2 

B (4800-5600 K) 

2.4±0.1 

4.9±0.1 

6.8±0.1 

B (5600-6400 K) 

3.5±0.2 

7.8±0.3 

7.8±0.2 

B (6400-7200 K) 

5.1±0.6 

10.3±1.1 

9.2±1.0 

Parameters defining the scaling function 5„sini 

(Eq. 4) 

a 

— 

0.047±0.005 

— 

C(kms“*) 

— 

15.8 

— 


Average value of [} for the mix of stars analysed in this paper 


p 0.63 

Variation of p with template temperature (section 4.4) 


j8 (3200-4000 K) 

— 

0.51 ±0.02 

— 

yS (4000-4800 K) 

— 

0.60±0.04 

— 

yS (4800-5600 K) 

— 

0.69±0.08 

— 

yS (5600-6400 K) 

— 

0.80±0.06 

— 

yS (6400-7200 K) 

— 

1.22±0.15 

— 


* Where C is calculated from Eq. 4 assuming a limb darkening 
coelficient u = 0.6 (Claret Diaz-Cordoves & Gimenez 1995). 


This is not unexpected since measurements of v sin i should be 
much less effected by binarity. A Student’s t-distribution with 
V = 2 fits either the short- or long-term repeat CDFs equally 
well. 

There are too few stars in our sample with v sin i > 5 km s“^ 
for a detailed investigation of how p might vary with age, log g 
or metallicity subsets. 

5. Measurement uncertainties using different 
instrumentai configurations 

So far the analyses have been restricted to observations with 
the HR15N order-sorting filter. In this section we consider 
how these results can be extended to the other GES observa¬ 
tional setups. We used all of the “GES_MW” (GES Milky Way 
Programme) fields, consisting of more than 20,000 RV measure¬ 
ments from individual spectra taken with the HR 10 and HR21 
filters. Unfortunately there are too few measurements through 
these filters with vsin/ > 5 km s ' to constrain the vsin/ preci¬ 
sion in the same way that was done for HR15n observations. 

These precision of the HR 10 and HR21 RV measurements 
were compared with those predicted using the simple model de¬ 
scribed in Appendix A. Eor this comparison it is assumed that: 

- The uncertainty in RV precision scales as Srv (See Sect. 3.1) 

- Parameter C characterising the dependence of RV un¬ 
certainty on vsin/ depends on the spectral resolution as 
0.895c/Ri (see Eq. A4). 

- Parameter A that determines the difference in RV precision 
for short- and long-term repeats corresponds to a displace¬ 


ment of the spectrum on the detector in the dispersion direc¬ 
tion, measured in pixels, rather than a fixed velocity differ¬ 
ence. The number of physical CCD pixels contributing each 
spectrum along the dispersion direction is 4096, so = 
4096/L4/cA/l where Ad is the wavelength range of the filter 
(see Table 3). In the case of filter HR15, A - 0.25 km s * 
corresponds to - 0.061 pixels, which we will assume is 
the same for the other filters. 

Only constant B has to be found, and this can be done using the 
distribution of Erv found from short-term repeats. This allows 
the measurement precision of RV for a given filter to be esti¬ 
mated even when there are no long-term repeat measurements 
to make an independent empirical analysis. In the case of filters 
HR 10 and HR21 it turns out that there are enough long-term re¬ 
peat measurements, albeit over a restricted range of v sin / values, 
to test this hypothesis. 

Figure 10 shows an analysis for all field stars that were ob¬ 
served in the GES_MW fields. Figures 10a and d show the vari¬ 
ation of the standard deviation of Erv x S fN with v sin i. Data 
with large v sin / values are few; therefore the error bars become 
large with increasing v sin i. Even so, the curve corresponding 
to the value of B evaluated for the full data set using the value 
of C predicted from Eq. A4 is consistent with the empirically 
measured uncertainties. 

Figures 10b and e show how the standard deviation of 
(Erv)(S /N)/(l + ([vsin/]/C)^)^^^ (estimated using the MAD) 
varies with S/N. For S/N < 100 both plots show reasonable 
agreement (within 10 percent) between the measured data and 
the line showing a single value of B evaluated for the full data 
set using theroretical C. Agreement is less good for data with 
S/N > 100. However, any inaccuracy here will have little ef¬ 
fect on the estimated uncertainty in RV for the majority of stars 
which are slow rotators since, at high values of S /N, the uncer¬ 
tainty of these stars is dominated by the constant term, A in the 
expresssion for S rv (see eqns. 1 and 2). 

The mean values of B are 2.3 km s"* for filter HR 10 and 
7.1 km s ' for filter HR21, compared with 5.0km s ' forHR15N 
i.e for spectra with similar vsin / and S/N RVs estimated from 
spectra taken with HRIO are more precise. The temperature de¬ 
pendence, illustrated in Fig. 5, is also different in detail. Figures 
6a and b show the variation of B with logg and [Fe/H]. The 
trends are similar to the variation found for order sorting filter 
HR15N i.e B is almost independent of gravity and changes only 
slowly with metallicity. 

Parameter A was determined in two ways. First, it was deter¬ 
mined using the measured data for the relatively small sample of 
long-term repeats as described in Sect. 3.1. This gave values of 
A = 0.18 + 0.02km s-' for filter HRIO and 0.28 + 0.02km s"' 
for filter HR21. These compare with the predicted values of 
0.22 km s“' and 0.26 km s“' inferred by scaling the dpix value 
determined for filter HR15N by the ratio of their pixel sizes in 
km s"'. 

Figures 10c and lOf shows the CDFs of the normalised un¬ 
certainty for short- and long-term repeats using filters HRIO and 
HR21 respectively. In each case S rv is evaluated using the ap¬ 
propriate theoretical values of A and C and the mean empirical 
value of B determined for each filter, and these are inserted into 
Eqns. 1 (for short-term repeats) or 2 (for long-term repeats). Also 
shown is the CDF of a Student’s t-distribution with v = 6 which, 
as for the HR15N data, appears to be an excellent representa¬ 
tion of the distribution due to short-term repeats. The data are 
sparse for long-term repeats, but the distributions appear to have 
more extended tails, consistent with the idea that they contain 
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Fig. 10. Analysis of the empirical uncertainty for short term repeat observations of RV using filter HR 10 and HR21. Plots a and d 
show the variation of Erv with vsini. Red lines show the curves predicted using the model value of C (see Eq. A4). The dashed 
lines shows the curves predicted using values of B and C fitted to the binned data. Plots b and e show the variation in the estimated 
value of parameter B in Eq. 1 with S /N. Red lines show values of B fitted to the full data set for each filter. Blue crosses show the 
estimated values B as a function of S /N corrected for the measured variation of B with Tefr (see Table 3). Plots c and f shows the 
normalised uncertainty for short term repeats and long term repeats. The black curve shows the CDE for short-term repeats. The red 
line shows the CDE for long-term repeats. Blue diamonds show a Student’s t-distribution with v = 6, which matches the distribution 
for short-term repeats well. 


RV shifts due to binary systems. We do not have sufficient data to 
test whether the uncertainty CDEs for long-term repeats within a 
night are similar to those for short-term repeats, but we assume 
that, like the HR15N data, this will the case for data taken with 
HRlOand HR21. 


6. Discussion and Summary 

We have shown that the normalisation functions given in Eqs. 
(1) -(4) are reasonable descriptions of how uncertainties in RV 
and V sin i scale with S/N and v sin i. The recommended average 
parameters of A, B, C defining the scaling function for RV are 
given in Table 3 for observations performed with the three main 
instrumental configurations used for GIRAEEE observations in 
the GES. Average values of a, yS, and C that define the scaling 
function for vsini are also given for filter HR15N. The uncer¬ 
tainties given by Eqns. 2 and 4 are not normally distributed; they 
have more extended tails. The uncertainty distribution for a given 
observation of RV is better represented by the value of Eq. 2 
multiplied by a normalised Student’s t-distribution with v = 6, 


whilst for V sin i the uncertainty distribution can be approximated 
by Eq. 4 multiplied by a normalised Student’s t-distribution with 
v = 2. 

Equations 2 and 4 decouple the influences of spectral type 
and the spectrograph; A, C and a are properties of the instru¬ 
mental setup, whilst B and (3 depend on the type of star observed. 
The dependence on gravity, age and metallicity, over the range 
-1 <[Ee/H]< 1, is weak; but the temperature dependence be¬ 
comes strong for T^ff > 5200 K, such that B and /3 increase with 
Teff and the precision worsens. This is presumably the result of a 
decreasing number of strong, narrow lines in the spectra of hot¬ 
ter stars. The temperature dependent B and fi values are listed in 
Table 3 and should be used in conjunction with the mean values 
of A, C and a. There are insufficient observations of stars with 
vsini > 5 km s ' using order-sorting filters HRIO and HR21, 
so we cannot estimate /? for such observations. It should also be 
noted that for reasons of sample size, the calibration of B and (3 
is limited to 3200 < < 7200 K. 

Parameter A is between 0.22 and 0.26 kms“', dependent on 
instrumental setup, and represents the best precision with which 
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RV can be obtained from an individual GES spectrum with low 
rotational broadening and large S IN. The origin of this term is 
unclear; it partly arises from uncertainties in the wavelength cal¬ 
ibration and the application of calibration offsets from the “sim- 
cal” fibres or sky emission lines. However, various tests have 
shown these cannot be entirely responsible and we suspect there 
are additional contributions that may be associated with move¬ 
ment of the fibres at the spectrograph slit assembly or target mis- 
centering in the fibres combined with imperfect signal scram¬ 
bling. 

The analyses we present were derived from results in the 
iDR2 GES data release and the coefhcients in Table 3 are appli¬ 
cable to those data and also to the more recent iDR3 update that 
used the same pipeline analysis. The exact use of these results 
depends on the purpose of any particular investigation. To esti¬ 
mate what approximates to a particular confidence interval on a 
RV (v sin i) value, the following procedure is recommended; 

1. Use the instrumental setup and an estimated stellar temper¬ 
ature (preferably from the GES analysis) to choose the ap¬ 
propriate values of B (JI) and C from Table 3; calculate S rv.q 
(•^vsim o) from Eq. 1 (Eq. 3) using the measured vsin/ and 
S/N. ’ 

2. Choose the A (a) value appropriate for the instrumental setup 
from Table 3 and calculate Srv (^vsim) from Eq. 2 (Eq. 4). 

3. If combining results from repeated observations, these 

should be weighted using short-term re¬ 

peats (i.e. without the inclusion of the A (a) term), or 5^^ 
(5^^jn,) for long-term repeats. 

4. Eor accurate modelling of RV data one should use Srv 
(S ysini) multiplied by a Student’s t-distribution with v = 6 
(v = 2) as a probability distribution for the uncertainty. More 
crudely, a confidence interval can be estimated by multi¬ 
plying SRV (Svsini) by the appropriate percentile point of a 
Student’s t-distribution with v = 6 (v = 2). Eor example, to 
estimate a 68.3 per cent error bar, multiply by 1.09 (1.32), or 
for a 95.4 per cent error bar multiply by 2.51 (4.50). 

Note that whilst the 68.3 per cent confidence intervals are 
quite close to the value expected for a normal distribution with a 
standard deviation of S rv (S vsini)^ the 95.4 per cent confidence 
intervals are significantly larger due to the broader tails of the 
Student’s t-distributions. We do not recommend extrapolating 
these estimates to even larger confidence intervals since we have 
few data with which to reliably constrain the distribution at these 
values. It seems likely that at the conclusion of GES there will 
be sufficient data (roughly 5 times as much) to significantly im¬ 
prove this situation. A larger dataset will also allow us to study 
how V sin i precision varies with and log g and between dif¬ 
fering observational setups. 
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Appendix A: Variation of measurement precision 
with radiai and projected rotation veiocities 


We consider below how the measurement precision of RV and 
V sin i scale with S /N and v sin i for short-term repeats where 
there are no changes in setup or wavelength calibration between 
observations. We make the simplifying assumption that the pre¬ 
cision in RV scales as Erv oc I(SIN) where W, is the 
EWHM of a Gaussian profile representing the characteristic ab¬ 
sorption line profile of the measured spectrum. This approximate 
relation can be deduced from the results of Landman, Roussel- 
Dupre and Tanigawa (1982). These authors showed that, for the 
ideal case of a Gaussian line profile of amplitude a, mean value 
m, and standard deviation s, sampled using binned data with a 
uniform Gaussian noise of rms amplitude e per bin, the statis¬ 
tical uncertainties in the estimated values of m and s are given 
by; 



(A.l) 


where is the uniform bin width and << s. 

In the present case a varies with equivalent width, (EW), of 
the characteristic absorption line as a = (EW)hl yfhi s, where h 
is the amplitude of the contiuum. If the depth of the absortion 
line, a, is small compared to the continuum, h then measurement 
uncertainty is e ^ hl(SIN). Substituting these values in equation 
Al using the relation W - V8 In 2 s gives; 


Ax1/2\JV3/2 /Ax1/2\^3/2 

- -and cTw oc - - 

EW j SIN \EW j SIN 


(A.2) 


A. 1. Effect of V sin i on FWHM of the absorption line 

Eor a slowly rotating star, assuming that any sources of broaden¬ 
ing other than rotation are much smaller than the intrinsic spec¬ 
trograph resolution, the EWHM of an individual absorption line 
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is Wo = where A is the mean wavelength and the resolv¬ 
ing power of the spectrograph. For fast rotating stars the width 
of the spectral lines is increased by rotational broadening. Gray 
(1984) gives the rotational broadening kernel as 


K(A)^ 


1 

A 


2(1 - U)I7T 

(I-m/3) 



m /2 

(l-«/3) 





where A = /l[vsin/]/c, A is wavelength (over the range -A < 
/I < A) and u is the limb darkening coefficient. 

Convolving a spectrum with this kernel increases the FWHM 
of individual lines approximately as W ^ sj^o + 
where Arms is the rms of the broadening kernel (A^^^ - 

J A^KdA). Evaluating jrms from Eq. A3 gives; 


W = Wo 


1 -H 


V sin i 
C 


2\ 


1/2 


(A.4) 


where C = 


/ i-m/3 
U-7h/15/ 


c 



A.2. Scaling of uncertainty in RV and v sin i 

To determine how the uncertainty in radial velocity, Erv scales 
with S/N and v sin i we assume Erv cTm- For a given spectra 
Ax and EW are independent of W and SIN so that (from eqns. 
A2 and A4) Erv scales with v sin i and S /N as , 


^RV,Q 


= B 


(l -H([vsin/]/C)2) 
S/N 


3/4 


(A.5) 


where B is an empirically determined constant and C depends 
on Rrt and u. A value of m = 0.6 is used in this paper (Claret 
Diaz-Cordoves & Gimenez 1995) giving C = 0.895c//?i. 

The uncertainty in the estimated value of v sin i is determined 
from the uncertainty in the estimated absorption line width, crw 
(eqns. A2 and A4) as; 


WV sin i 


CWo-w 

W2 - 1 


(A.6) 


Using this expession the uncertainty in the normalised value of 
Vsin i, ( oc cTj, sin,7[v sin /]) scales with vsin i and S/N as; 


5ysin/,0 “ 


(l+([vsin/]/C)2f" 
(5/A)([vsin/]/C)2 ’ 


(A.7) 


where p is an empirically determined constant. 
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